A  STATISTICAL  MODEL  OF  ATMOSPHERIC 
TEMPERATURE  SIGNALS 


Edward  Marvin  Kline 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


TH 


A  STATISTICAL  MODEL  OF 
ATMOSPHERIC  TEMPERATURE  SIGNALS 


by 


Edward  Marvin  Kline 


March  1972 


AppJiov&d  ioh.  public  h.oXzjte<L',    (LUthAhuutLovi  unlAjnittd. 


A  Statistical  Model  of 
Atmospheric  Temperature  Signals 


by 


Edward  Marvin  Kline 
Lieutenant,  United  States  Navy 
BoS.,  University  of  Idaho,  1967 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OCEANOGRAPHY 

from  the 
waval  POSTGRADUATE  SCHOOL 


ABSTRACT 

The  "ramp"  is  an  often  observed  feature  in  temperature 
fluctuations  during  unstable  atmospheric  conditions.   It  is 
characterized  by  a  gradual  increase  in  temperature  followed 
by  a  sudden  drop  to  an  ambient  level.   This  ramp  clearly 
distinguishes  temperature  signals  from  other  turbulence  signals 
such  as  velocity.   Three  different  ramp-type  atmospheric 
temperature  fluctuations  and  their  derivatives  are  constructed 
and  statistically  examined  for  the  parameters  skewness  and 
coefficient  of  excess.,   These  statistical  values  are  compared 
with  values  obtained  from  actual  signals.   The  "complex  ramp 
signal"  was  found  to  represent  these  quite  well.   Skewness  and 
coefficient  of  excess  for  this  signal  could  be  altered  almost 
independently  of  each  other.   The  complex  ramp  is  constructed 
by  inscribing  nearly  symmetrical  consecutive  triangles  within 
two  envelope  rays  of  constant  but  different  slopes  such  that 
the  two  envelope  rays  approximate  a  ramp  function. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION  8 

A.  OBJECTIVE   .  .  .  . 8 

B.  GENERAL  FEATURES  OF  TEMPERATURE  SIGNALS   ...  9 

II.  SIGNALS 13 

A.  ATMOSPHERIC  TEMPERATURE  FLUCTUATIONS  13 

B.  MODELS  OF  ATMOSPHERIC  TEMPERATURE  FLUCTUATIONS  13 

III.  METHOD  OF  EXAMINATION 16 

A.  SIGNAL 16 

B.  DIFFERENTIATED  SIGNALS 18 

C.  STATISTICS  APPLIED  19 

IV.  RESULTS 26 

V.  CONCLUSIONS 49 

VI.  RECOMMENDATIONS  FOR  FURTHER  STUDY 50 

COMPUTER  PROGRAM    53 

BIBLIOGRAPHY   59 

INITIAL  DISTRIBUTION  LIST   61 

FORM  DD  1473 62 


FIGURE  1: 
FIGURE  2: 
FIGURE  3: 
FIGURE  4: 
FIGURE  5: 
FIGURE  6: 
FIGURE  7: 
FIGURE  8: 
FIGURE  9: 
FIGURE  10: 
FIGURE  LI: 
FIGURE  12: 
FIGURE  13: 
FIGURE  14: 
FIGURE  15: 


LIST  OF  FIGURES 

■ 
Atmospheric  Temperature  and  Differentiated 
Temperature  Signal.  ..„ ]_]_ 

Schematic  of  Convective  Plume  According  to 

Gill  (1971) 12 

Growth  of  Ramps  as  Atmosphere  Changes  from 
Stable  to  Unstable 15 

Standard  Ramp  Model  and  Differentiated 

Standard  Ramp  Model  ....„  22 

Asymmetrical  Ramp  Model  and  Differentiated 
Asymmetrical  Ramp  Model  23 

Complex  Ramp  Model  and  Differentiated 

Complex  Ramp  Model  „ 24 

Variation  of  Standard  Ramp  Signal  with 

Period  ...„ 28 

Variation  of  Standard  Ramp  Signal  with 

Ramp  Width   „ „ 31 

Variation  of  Standard  Ramp  Signal  with 

Ramp  Slope 0 31 

Variation  of  Asymmetrical  Ramp  Signal  with 
Period 33 

Variation  of  Asymmetrical  Ramp  Signal  with 

Ramp  Width «, 34 

Variation  of  Asymmetrical  Ramp  Signal  with 

Back  Slope  .  , 34 

Variation  of  Asymmetrical  Ramp  Signal  with 

Back  Slope  Width 38 

Variation  of  Complex  Ramp  Signal  with 

Period „   »  .  . 40 

Variation  of  Complex  Ramp  Signal  with 

Upper  Basic  Slope 42 


FIGURE  16:   Variation  of  Complex  Ramp  Signal  with 

Lower  Basic  Slope  4Z4. 

FIGURE  17:   Variation  of  Complex  Ramp  Signal  with 

Ramp  Width  ..„ •  46 


LIST  OF  TABLES 


TABLE  I. 
TABLE  II. 
TABLE  III.   - 
TABLE  IV. 
TABLE  V. 
TABLE  VI. 
TABLE  VII.   - 
TABLE  VIII.  - 
TABLE  IX. 
TABLE  Xo 
TABLE  XI. 
TABLE  XII.   - 
TABLE  XIII.  - 
TABLE  XIV.   - 


SKEWNESS  AND  COEFFICIENT  OF  EXCESS  OF 
DIFFERENTIATED  TEMPERATURE  SIGNALS  ....   25 

VARIATION  OF  DIFFERENTIATED  RAMP  SIGNAL 

WITH  PERIOD 27 

VARIATION  OF  DIFFERENTIATED  RAMP  SIGNAL 

WITH  RAMP  SLOPE 29 

VARIATION  OF  DIFFERENTIATED  RAMP  SIGNAL 

WITH  RAMP  WIDTH 30 

VARIATION  OF  DIFFERENTIATED  ASYMMETRICAL 

RAMP  SIGNAL  WITH  PERIOD 32 

VARIATION  OF  DIFFERENTIATED  ASYMMETRICAL 

RAMP  SIGNAL  WITH  BACK  SLOPE 34 

VARIATION  OF  DIFFERENTIATED  ASYMMETRICAL 

RAMP  SIGNAL  WITH  RAMP  WIDTH 35 

VARIATION  OF  DIFFERENTIATED  ASYMMETRICAL 

RAMP  SIGNAL  WITH  BACK  SLOPE  WIDTH   ....   37 

VARIATION  OF  DIFFERENTIATED  COMPLEX  RAMP 
SIGNAL  WITH  PERIOD 39 

VARIATION  OF  DIFFERENTIATED  COMPLEX  RAMP 
SIGNAL  WITH  UPPER  BASIC  SLOPE   41 

VARIATION  OF  DIFFERENTIATED  COMPLEX  RAMP 
SIGNAL  WITH  LOWER  BASIC  SLOPE   43 

VARIATION  OF  DIFFERENTIATED  COMPLEX  RAMP 
SIGNAL  WITH  RAMP  WIDTH 45 

VARIATION  OF  DIFFERENTIATED  COMPLEX  RAMP 
SIGNAL  WITH  PERIOD  AND  RAMP  WIDTH   ....   47 

VARIATION  OF  A  DIFFERENTIATED  VELOCITY 

SIGNAL  WITH  RAMP  WIDTH 52 


ACKNOWLEDGEMENT 

The  author  would  like  to  thank  Professor  N.  E.  J.  Boston, 
whose  appreciated  advice  and  encouragement  greatly  contributed 
to  this  thesis.   The  author  is  also  grateful  to  Professor 
K.  L.  Davidson  for  his  willing  assistance.   Miss  S.  D.  Raney 
of  the  W.  R.  Church  Computer  Center  provided  much  needed  pro- 
gramming advice  and  Professor  A.  L.  Schoenstadt  provided 
mathematical  assistance.   Finally,  the  author  would  like  to 
thank  his  wife,  Linda,  who  served  as  part-time  typist  and 
proof  reader o 


I.   INTRODUCTION 

A.   OBJECTIVE 

Atmospheric  and  oceanic  turbulence  contributes  greatly  to 
interactions  at  the  air /sea  interface.  Mathematical  models 
have  been  developed  which  attempt  to  describe  this  turbulent 
boundary  layer  where  atmosphere  and  ocean  meet.   As  sensing 
systems  become  more  sophisticated  and  more  turbulence  data 
becomes  available,  existing  models  are  accordingly  modified. 
An  important  modification  was  made  by  Kolmogorov  (1962)  and 
others  to  account  for  rather  large  fluctuations  observed  in 
the  energy  dissipation  rate0   Such  fluctuations  were  observed 
primarily  in  differentiated  turbulence  signals  which  emphasized 
high  frequency  components.   Kolmogorov  revised  his  theory  in 
order  to  incorporate  how  the  velocity  derivative's  skewness 
(s)  and  coefficient  of  excess  (c.e.)  increased  with  the 
Reynold's  number  (Re).   In  the  energy  containing  region  the 
parameters  s  and  c.e.  are  zero  since  the  turbulence  fluctuations 
seem  to  have  a  Gaussian  distribution.  At  small  scales  (high 
frequencies)  they  are  decidedly  not  zero.   The  skewness  and 
coefficient  of  excess  of  the  derivative  of  the  temperature 
fluctuations  are  not  zero  either.   One  might  then  ask,  "Knowing 
s  and  c.e.  from  real  data,  is  it  possible  to  construct  signals 
that  resemble  observed  temperature  signals?"   Such  was  the. 
objective  of  this  thesis. 
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B.   GENERAL  FEATURES  OF  TEMPERATURE  SIGNALS 

The  "ramp"  is  an  often  observed  feature  in  temperature 
fluctuations  during  unstable  atmospheric  conditions.   It  is 
characterized  by  a  gradual  increase  in  temperature  followed 
by  a  sudden  drop  to  an  ambient  level.   The  signal  rarely  falls 
below  this  ambient  level  (Figure  1).   This  sudden  drop  clearly 
distinguishes  a  temperature  signal  from  other  turbulence 
signals  such  as  velocity  signals  whose  main  perturbations 
from  the  ambient  level  are  usually  gradual  rises  with 
accompanying  gradual  descents.   It  seems  natural  to  conclude 
that  these  observable  non-Gaussian  type  features  in  temperature 
signals  probably  contribute  to  the  measureable  skewness  and 
coefficient  of  excess. 

Electronic  pattern  recognition  studies  of  atmospheric 
temperature  fluctuations  have  been  made  by  Sokol  (1971). 
Although  his  study  was  not  exhaustive,  he  achieved  over  70% 
detection  success.   Another  method  of  describing  these  signals 
is  statistically  through  skewness  and  coefficient  of  excess. 
Since  these  statistics  are  zero  for  Gaussian  noise,  their 
values  would  represent  only  the  non-Gaussian  character  of  the 
signal.   Thus,  it  may  be  possible  in  the  future  to  use  the 
frequency  and  structure  of  the  ramps  as  a  measure  of  some 
meteorological  conditions. 

Kaimal  and  Businger  (1970)  and  Gill  (1971)  have  demonstrated, 
over  land  at  least,  that  the  ramps  are  manifestations  of 
convective  plumes  present  near  the  air/earth  interface.   A 
schematic  of  the  plume  and  its  relation  to  the  ramp,  by 
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Gill  (1971),  are  shown  in  Figure  2.   Taylor  (1958)  and  Gill 
(1971)  have  shown  that  these  ramps  can  be  correlated  to  other 
ramps  appearing  simultaneously  at  different  heights  but  with 
a  slight  offset  owing  to  a  tilt  in  the  plume;  presumably  due 
to  wind  shear.   Boston  (1971)  indicated  that  the  ramps  in  the 
temperature  signal  often  correlate  with  perturbations  in 
velocity  signals,   Safley  (1972)  has  related  them  to  water 
waves  and  velocity  signals. 

Plumes  form  a  significant  feature  in  the  atmospheric 
boundary  layer.   Although  they  are  present  only  a  small 
fraction  of  the  time  they  may  account  for  the  majority  of  heat, 
moisture,  and  momentum  exchange.   Safley  (1972)  found  that 
the  plumes  accounted  for  10%  of  the  signal  but  32%  of  the 
heat  flux.   Furthermore,  they  cause  random  fluctuations  in 
electromagnetic  waves  over  line-of-sight  paths.   Of  recent 
interest  is  the  effect  they  have  on  laser  beams.   When  beamed 
over  water  they  exhibit  a  random  walking  at  the  target.   This 
is  of  considerable  importance  if  a  laser  beam  is  to  be  used 
in  a  variety  of  practical  applications ,  such  as  fire  control 
systems . 
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II.   SIGNALS 

A.  ATMOSPHERIC  TEMPERATURE  FLUCTUATIONS 

When  the  atmosphere  is  stable  the  temperature  fluctuations 
tend  to  be  small  because  the  vertical  turbulent  velocity  com- 
ponent tends  to  be  damped  out  by  the  vertical  temperature 
gradient.   When  the  atmosphere  is  neutral  the  fluctuations 
are  also  small  but  tend  to  have  both  positive  and  negative 
excursions  due  to  the  lack  of  restoring  buoyant  forces. 

Unstable  conditions,  on  the  other  hand,  exhibit  the 
characteristic  ramp  features  in  the  temperature  trace  (Figure  3). 
This  is  the  condition  that  usually  exists  during  the  daytime 
in  the  lower  atmosphere  because  the  earth  is  heated  and  warms 
from  below  causing  upward  sensible  heat  transfer. 

B.  MODELS  OF  ATMOSPHERIC  TEMPERATURE  FLUCTUATIONS 
Atmospheric  turbulence  for  stable  and  neutral  conditions 

is  difficult  to  model  since  it  is  a  broad  band  stochastic 
process.   However,  because  of  the  characteristic  ramp  shape 
of  temperature  fluctuations  in  unstable  conditions,  some 
attractive  possibilities  for  models  present  themselves.   The 
model  for  the  temperature  signal  of  unstable  conditions  will 
consist  of  a  basically  straight  line  ramp  signal  to  which 
variations  will  be  added 0 

Rather  than  examining  spectra  of  the  modeled  atmospheric 
temperature  fluctuations,  statistics  of  the  derivative  of  the 
model  will  be  examined.   The  differentiated  signal  emphasizes 
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the  high  frequency  (high  wave  number)  region  of  the  signal 
which  is  the  region  of  most  interest  in  the  modeling  process 
since  measurements  in  the  last  decade  clearly  indicate  that 
the  high  wave  number  region  is  decidedly  non-Gaussian.  Also, 
visual  inspection  of  ramps  indicates  that  high  frequency 
fluctuations  may  be  important  in  ramp  structures. 
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III.   METHOD  OF  EXAMINATION 

A.   SIGNAL 

The  signal  modeling  was  done  on  an  IBM  360/67  digital 
computer  at  the  W.  R.  Church  Computer  Center,  Naval  Postgraduate 
School,  Monterey,  California.   The  basic  ingredient  in  each 
model  is  the  individual  ramp  which  will  be  repeated  at  a 
specified  interval.   This  collection  of  repetitive  ramps  will 
comprise  the  model  signal.   Each  signal  will  be  stored  in  a 
10,000  X  1  vector  array  which  represents  the  time  axis  of  the 
signal.   The  size  of  the  array  may  be  varied,  however,  since 
the  models  are  periodic.   The  distance  between  each  array 
entry  is  constant  and  equal  to  the  distance  between  samples  of 
a  digital  recorder  or  converter  (which  is  a  function  of  the 
original  sampling  rate). 

The  models  are  dimensionless  although  voltage  is  implied 
and  terms  such  as  period  and  width  must  be  redefined  as  they 
refer  to  digital  words.  Width  refers  to  consecutive  digital 
words  or  array  entries  and  is  used  to  describe  the  length  of 
a  ramp  or  a  ramp  component.  Ramp  period  is  the  number  of 
consecutive  digital  words  between  the  start  of  a  ramp  and  the 
start  of  the  next  succeeding  ramp. 

Three  ramp  modifications  of  increasing  complexity  will  be 
examined  as  possible  models  for  the  temperature  signal.   The 
first  is  a  "standard  ramp  signal"  whose  ramps  are  constructed 
of  a  gradual  rise  ending  in  a  vertical  drop  (Figure  4)  as 
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described  by  Gibson  e_t  ai  (1971).   The  variable  parameters  of 
this  signal  are  (1)   the  period  of  the  ramp,  (2)   the  slope 
of  the  ramp,  and  (3)   the  width  of  the  ramp. 

The  second  modification  is  an  "asymmetrical  ramp"  with  a 
gradual  rising  front  slope  followed  by  a  short  steep  back 
slope  (Figure  5).   The  variables  in  this  signal  are  (1)   ramp 
period,  (2)   the  back  slope  of  the  ramp,  (3)   ramp  width,  and 
(4)   the  width  of  the  back  slope.   The  front  slope  value  and 
amplitude  need  not  be  included  since  they  are  functions  of 
the  other  variables.   Specifically, 

front  slope  =  ((back  slope  width)  X  (back  slope))/ 

((ramp  width)  -  (back  slope  width)),  amplitude  = 

(back  slope)  (back  slope  width). 

The  third  modification,  a  "complex  ramp",  is  constructed 
by  inscribing  nearly  symmetrical  consecutive  triangles  (teepees) 
within  two  envelope  rays  of  constant  but  different  slopes  such 
that  the  two  rays  approximate  a  ramp  function  (Figure  6). 
Each  succeeding  slope  of  the  inscribed  teepees  will  be  greater 
in  absolute  value  and  alternate  in  sign.   This  type  of  model 
is  chosen  because  close  inspection  of  temperature  strip  charts 
shows  that  large  fluctuations  of  increasing  amplitude  is  a 
commonly  observed  ramp  feature.   The  variables  for  this 
signal  are  (1)   ramp  period,  (2)   slope  of  the  upper  defining 
ray,  (3)   slope  of  the  lower  defining  ray,  and  (4)   the  width 
of  the  ramp.   The  period  of  the  internal  teepees  is  not  listed 
as  a  variable  because,  as  will  be  explained  later,  the  model 
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is  dimensionless  and,  for  example,  halving  the  period  or 
doubling  the  frequency  of  the  internal  teepees  can  be  accom- 
plished by  simultaneously  doubling  the  ramp  period  and  ramp 
width, 

B.   DIFFERENTIATED  SIGNAL 

As  previously  mentioned  the  differentiated  signal  will  be 
examined  since  it  emphasizes  the  higher  frequencies  that 
contribute  to  the  non-Gaussian  nature  of  the  signals.   There- 
fore, now  that  the  direct  signal  has  been  modeled  the 
differentiated  signal  must  be  computed.   The  differentiated 
standard  ramp  signal  is  merely  a  box  function  whose  amplitude 
is  equal  to  the  slope  of  the  ramp  (Figure  4). 

The  differentiated  asymmetrical  ramp  signal  is  an  alternating 
asymmetrical  box  function.   Each  differentiated  pulse  consists 
of  a  broad  positive  box  whose  amplitude  is  equal  to  the  front 
slope  of  the  ramp  followed  by  a  short,  highly  negative  box 
whose  amplitude  is  equal  to  the  back  slope  of  the  ramp  (Figure  5). 

The  signal  model  for  the  derivative  of  the  complex  ramp 

(Figure  6)  is  constructed  by  letting 

dy  =  us lope  (1),  the  slope  of  the  uppermost  limiting  ray  and, 
dx 

dy_'=  Is  lope  (2),  the  slope  of  the  lowermost  limiting  ray. 
dx1 

For  representation  on  a  digital  computer  the  implied  distance 
between  successive  array  entries  is  a  mesh  length.   Then  within 
the  ramps  each  succeeding  array  entry  represents  successive 
slopes  of  the  inscribed  teepees.   However,  for  purposes  of  com- 
puting these  slopes  the  mesh  length  is  offset  by  half  its 
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length  and  the  difference  between  successive  intersections 
of  the  upper  slope  and  the  lower  slope  is  computed.   From 
equations  (1)  and  (2), 

dy  =  (dx)(uslope) 

dy'  =  (dx') (lslope). 
Let  n  =  dx  ,  where  n  is  the  number  of  mesh  lengths  along  the 
time  axis  of  the  signal  measured  from  the  start  of  the  ramp. 
Hence, 

dy  =  n(uslope) 

dy'n  =  n(lslope). 
Therefore,  assuming  the  first  slope  is  positive: 

x(n)   =  (n+1)  (uslope)  -  (n)  (lslope) 

x(n+l)=  (n+2)  (lslope)  -  (n+1)  (uslope) 

x(n+2)=  (n+3)  (uslope)  -  (n+2)  (lslope) 

x(n+3)=  (n+4)  (lslope)  -  (n+3)  (uslope) 

x(n+4)=  (n+5)  (uslope)  -  (n+4)  (uslope) 


Note  that  the  first  and  last  values  may  have  to  be  computed 
separately  since  the  signal  must  be  continuous  with  the  ambient 
level  at  the  beginning  and  end  of  the  ramp.   A  sample  program 

with  this  representation  is  given  in  a  separate  section. 

\ 

C.   STATISTICS  APPLIED 

One  can  usually  discern  a  temperature  signal  from  a  velocity 
signal  by  the  ramp.   However,  if  the  spectra  of  these  signals 
are  examined  in  the  range  of  about  1  Hz  to  1000  Hz  they  are 
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remarkably  similar,  although  their  statistics  such  as  skewness 
and  coefficient  of  excess  are  quite  different.   As  a  result, 
recent  turbulence  literature  has  used  these  two  parameters 
extensively  to  describe  turbulence  signals.   Therefore,  skew- 
ness and  coefficient  of  excess  will  be  used  as  modeling 
parameters.   The  skewness  and  c.e.  of  the  models  will  be  com- 
pared to  those  of  actual  signals. 

Skewness  is  a  measure  of  the  amount  a  signal  differs  from 
its  mean  value.   The  skewness  of  a  Gaussian  signal  is  zero. 
A  distribution  that  is  skewed  to  the  right  has  a  positive 
skewness  and  a  distribution  that  is  skewed  to  the  left  has 
a  negative  skewness.   A  symmetrical  signal  will  have  zero  skew- 
ness and  a  signal  with  a  high  skewness  value  will  have  pulses 
that  are  very  asymmetrical.   Coefficient  of  excess  is  a 
measure  of  the  peakedness  of  a  distribution.   A  distribution 
that  is  less  peaked  than  normal  distribution  has  c.e.  values 
less  than  zero,  and  vice  versa.   The  c.e.  of  a  signal  is 
highly  sensitive  to  infrequent  deviations  from  the  mean  signal 
value.   Therefore,  a  signal  with  ramps  spaced  relatively  far 
apart  will  have  a  higher  c.e.  value  than  a  signal  with  ramps 
spaced  close  together <,   These  higher  order  statistics  emphasize 
the  higher  frequencies  which  are  non-Gaussian.   However,  since 
they  are  exponential  functions,  they  require  a  longer  record 
to  achieve  stable  results.   Therefore,  the  use  of  statistics 
of  orders  higher  than  c.e.  is  somewhat  limited. 

Skewness  (s)  is  defined  as  : 

3    9   3/2 
s  =  e  /(e  )    ,  where  e  is  the  variation  of  the 

signal  from  its  mean  value. 
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Coefficient  of  excess  (c.e.)  is  defined  as: 

4  U        o    2 

c-e*  "  e   ,  -3.   Often  one  sees  e  /(e  )  ,  the  ratio 

(e2)z 
of  the  fourth  moment  to  the  second  moment,  used  as  a  statistic 
and  referred  to  as  "kurtosis,"  but  c.e.  will  be  used  here 
because  it  is  zero  for  a  Gaussian  signal. 

The  skewness  and  c.e0  for  differentiated  temperature 
signals  have  been  computed  by  several  people  (Table  I).   Al- 
though the  one  thing  the  statistics  seem  to  have  in  common  is 
their  variability,  representative  figures  can  be  chosen.   For 
temperature  signal  a  skewness  value  between  0.2  and  0.5  seems 
common  and  likewise  a  value  for  c.e.  is  about  30.   Compare 
these  values  to  representative  velocity  signal  values  of  -0.5 
and  15  for  s  and  c.e.  respectively. 
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TABLE  I. 

SKEWNESS  AND  COEFFICIENT  OF  EXCESS  OF 
DIFFERENTIATED  TEMPERATURE  SIGNALS 


INVESTIGATOR 
BOSTON  (1971) 


GIBSON  et.  al.  (1970) 


SKEWNESS 

-0, 

,013 

0, 

,195 

0. 

.140 

0. 

,002 

-0. 

,082 

0, 

,596 

0. 

,203 

0, 

.307 

0, 

,315 

0, 

,171 

0, 

.673 

0, 

,039 

0. 

,286 

0, 

,251 

0, 

.257 

0, 

,051 

0, 

,40 

0. 

,72 

0. 

.48 

0, 

.72 

C.E. 

20.3 
17.1 
21.4 
22.6 
13.6 
22.8 
12.2 
13.4 
32.0 
13.9 
30.7 
15.5 
18.1 
11.0 
12.3 
6.0 

40 
35 
34 
22 


GURVICH  (1967) 


15  - 
14001 
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IV.   RESULTS 

The  standard  ramp  signal  proved  to  be  a  poor  representation 
of  the  temperature  signal.   For  an  appropriate  c.e.  the  skew- 
ness  was  positive  but  much  too  large;  about  five.   This  agrees 
with  the  results  obtained  by  Gibson  et.al.(1971)  from  a 
similar  model.   The  asymmetrical  ramp  was  no  improvement.   With 
this  signal  a  large  positive  c.e.  was  accompanied  by  a  large 
negative  skewness.   This  signal  proved  to  be  of  some  value 
apart  from  its  intended  use  and  is  discussed  in  Section  VI. 

The  complex  ramp  signal  showed  more  favorable  results.   Not 
only  can  compatible  skewness  and  c.e.  be  found  but  the  signal 
can  be  adjusted  so  that  skewness  and  c.e.  vary  almost  indepen- 
dently.  A  more  detailed  analysis  of  the  individual  signals 
follows.   It  should  be  kept  in  mind  that  while  a  direct  signal 
model  is  the  final  product,  the  different  model  theories  were 
tested  by  the  statistics  of  the  differentiated  signal.   There- 
fore, the  following  results  apply  to  the  statistics  of  the 
differentiated  signal. 
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TABLE  II. 


VARIATION  OF  DIFFERENTIATED  RAMP  SIGNAL 
WITH  PERIOD  (FIGURE  7) 


PERIOD 

RAMP  SLOPE 

RAMP  WIDTH 

S. 

C.E. 

10 

0.408 

-1.83 

20 

\ 

1.45 

0.250 

30 

2.15 

2.65 

40 

2.0 

4.0 

2.63 

4.90 

80 

3.98 

13.8 

LOO 

4.70 

20.0 

200 

6.48 

40.0 

400 

8.49 

70.0 

Note  here  that  skewness  and  c.e.  both  increase  with 
increasing  period  but  they  never  simultaneously  achieve  the 
correct  values  required  by  the  models. 
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TABLE  III. 


VARIATION  OF  DIFFERENTIATED  RAMP  SIGNAL 
WITH  RAMP  SLOPE  (FIGURE  9) 


PERIOD 

RAMP  SLOPE 

RAMP  WIDTH 

S.' 

C.E. 

0.1 

2.67 

5.12 

0.5 

2.67 

5.12 

40 

1.0 

4.0 

2.67 

5.12 

5.0 

2.67 

5.12 

S.  and  c.e.  are  independent  of  ramp  slope. 
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TABLE  IV. 


VARIATION  OF  DIFFERENTIATED  RAMP  SIGNAL 
WITH  RAMP  WIDTH  (FIGURE  8) 


PERIOD 

RAMP   SLOPE 

RAMP  WIDTH 

S. 

C.E. 

2.0 

4.13 

15.1 

3.0 

3.23 

8.41 

40 

2.0 

4.0 

2.67 

5.11 

5.0 

2.27 

3.14 

6.0 

1.96 

1.85 

7.0 

1.71 

0.929 

Here  the  results  are  similar  to  those  obtained  with  a 
varying  period.  Although  values  of  s  and  c.e.  can  be  achieved 
with  the  standard  ramp  signal  that  are  respectively  commensurate 
with  measured  values  from  actual  signals,  the  values  of  s  and 
c.e.  are  not  simultaneously  compatible  with  measured  signals 
and  hence  this  signal  is  discarded  as  an  appropriate  model. 
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RAMP  WIDTH 


FIGURE  8.   Variation  of  Standard  Ramp  Signal 
with  Ramp  Width  (Period  =  40,  Ramp 
Slope  =  2.0) 
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FIGURE  9.  Variation  of  Standard  Ramp  Signal 
with  Ramp  Slope  (Period  =  40,  Ramp 
Width  =  4.0) 
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TABLE  V. 


VARIATION  OF  DIFFERENTIATED  ASYMMETRICAL 
RAMP  SIGNAL  WITH  PERIOD  (FIGURE  10) 


PERIOD 

BACK 
SLOPE 

FRONT 
SLOPE 

RAMP 
WIDTH 

BACK 
SLOPE  WIDTH 

S. 

C.E, 

10 

-3.77 

13.2 

30 

-4.64 

21.6 

50 

-6.00 

38.0 

100 

-1.5 

0.167 

10 

1 

-8.48 

78.9 

300 

-14.7 

243 

500 

-18.9 

403 

Although  c.e.  values  are  achieved  that  agree  with  actual 

signal  values,  the  skewness  values  are  not  only  too  large  in 

absolute  value  but  negative  instead  of  the  small  positive 
values  desired. 
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FIGURE  10. 


Variation  of  Asymmetrical  Ramp  Signal 
with  Period 

(Back  Slope  =  -1.5,  Ramp  Width  =  10, 
Back  Slope  Width  =  1) 
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TABLE  VI. 


VARIATION  OF  DIFFERENTIATED  ASYMMETRICAL  RAMP 
SIGNAL  WITH  BACK  SLOPE  (FIGURE  12) 


PERIOD 

BACK 

FRONT 

RAMP 

BACK  SLOPE 

S. 

C.E 

SLOPE 

SLOPE 

WIDTH 

WIDTH 

-0.10 

0.011 

-8.43 

78.1 

-0.50 

0.056 

-8.43 

78.1 

100 

-1.0 

0.10 

10 

1 

-8.43 

78.1 

-1.5 

0.144 

-8.43 

78.1 

-1.9 

0.189 

-8.43 

78.1 

Once  again  the  slope  of  the  signal  has  no  effect  on  the 
modeling  statistics. 
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TABLE  VII. 


VARIATION  OF  DIFFERENTIATED  ASYMMETRICAL  RAMP  SIGNAL 
WITH  RAMP  WIDTH  (FIGURE  II) 


PERIOD      BACK    FRONT    RAMP     BACK 

SLOPE     SLOPE   WIDTH  SLOPE  WIDTH 


S. 


C.E. 


100 


-1.5 


0.30 

6 

0.25 

7 

0.21 

8 

0.19 

9 

0.17 

10 

0.08 

20 

0.05 

30 

7.30 

67.0 

7.72 

70.8 

8.02 

73.8 

8.25 

76.2 

8.43 

78.1 

9.23 

87.3 

9.49 

90.4 

These  statistics  are  not  compatible  with  measured  values. 
The  skewness  is  too  large  in  absolute  value  and  negative  in 
sign. 
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RAMP  WIDTH 


FIGURE    11-  VARIATION  OF  ASYMMETRICAL  RAMP  SIGNAL  WITH  RAMP  WIDTH 
"(PERIOD  =  100,  BACK  SLOPE  =-1.5,  BACK  SLOPE  WIDTH  =  1) 


-2.0 


-1.0 


BACK  SLOPE 
FIGURE    12 -VARIATION  OF  ASYMMETRICAL  RAMP  SIGNAL  WITH  3ACK  SLOPE 
(PERIOD  -  100,   RAMP  WIDTH  -  10,    BACK  SLOPE  WIDTH  =  1) 
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TABLE  VIII 


VARIATION  OF  DIFFERENTIATED  ASYMMETRICAL  RAMP 
SIGNAL  WITH  BACK  SLOPE  WIDTH  (FIGURE  13) 


PERIOD 

BACK 

FRONT 

RAMP 

BACK  SLOPE 

S. 

C.E. 

SLOPE 

SLOPE 

WIDTH 

WIDTH 

0.167 

1 

-8.48 

78.9 

0.375 

2 

-4.77 

29.8 

LOO 

-1„5 

0.643 

10 

3 

-2.77 

14.8 

1.0 

4 

-1.30 

8.78 

1.5 

5 

0.0 

7.0 

Note  that  the  skewness  finally  becomes  non-negative  (zero) 
but  the  c.e.  is  too  low  and  the  shape  of  the  signal  is  that  of 
a  teepee,  not  a  rampa  Therefore,  due  to  the  negative  skewness 
inherent  in  this  signal,  it  is  discarded  as  a  possible  model. 
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VARIATION  OF  ASYMMETRICAL  RAMP  SIGNAL  WITH  BACK  SLOPE  WIDTH 
(PERIOD  -  100,  BACK  SLOPE  «  -1.5,  RAMP  WIDTH  -  10) 

FIGURE  13. 
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TABLE  IX. 


VARIATION  OF  DIFFERENTIATED  COMPLEX  RAMP  SIGNAL 
WITH  PERIOD  (FIGURE  14) 


LOWER  BASIC   RAMP    S.      C.E. 
SLOPE   WIDTH 


PERIOD 

UPPER  BASIC 

SLOPE 

40 

80 

100 

200 

300 

0.8 

400 

600 

800 

1000 

-0.4 


0.046 

-1.27 

0.064 

0.459 

0.072 

1.32 

0.102 

5.64 

40   0.124 

9.71 

0.145 

14.3 

0.175 

22.4 

0.201 

30.2 

0.229 

40.2 

This  signal  closely  approximates  actual  measured  values 
for  periods  in  the  range  of  800  to  1000.   This  appears  to  be 
an  excellent  model  for  typical  atmospheric  temperature  signals 
during  unstable  conditions.   Now  the  effect  of  the  remaining 
parameters  will  be  examined. 
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FIGURE  14 :  VARIATION  OF  COMPLEX  RAMP  SIGNAL  WITH  PERIOD 
(UPPER  BASIC  SLOPE  -0,8,  LOWER  BASIC  SLOPE  =  -0.U,  RAMP  WIDTH  «  10) 
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TABLE  X. 


VARIATION  OF  DIFFERENTIATED  COMPLEX  RAMP  SIGNAL 
WITH  UPPER  BASIC  SLOPE  (FIGURE  15) 


PERIOD    UPPER  BASIC 
SLOPE 

0.8 
1.2 
800  1.6 
2.0 
2.4 
50 


LOWER   BASIC 

RAMP 

S. 

C.E 

SLOPE 

WIDTH 

0.208 

30.6 

0.201 

30.2 

0.175 

30.2 

-0.4 

40 

0.153 

30.2 

0.134 

30.2 

0.119 

30.2 

0.008 

30.7 

The  coefficient  of  excess  remains  almost  constant  while 
the  skewness  ranges  over  values  not  unlike  those  encountered 
in  real  signals „ 
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TABLE  XI. 


VARIATION  OF  DIFFERENTIATED  COMPLEX  RAMP  SIGNAL 
WITH  LOWER  BASIC  SLOPE  (FIGURE  16) 


PERIOD     UPPER  BASIC     LOWER  BASIC   RAMP      S.       C.E. 
SLOPE  SLOPE     WIDTH 

0.6 

0.4 

0.2 
800  0.8  0.0      40 

-0.2 
-0.4 
-0.6 


-10.3 

227 

-  2.29 

58.2 

-  0.488 

34.6 

0.0 

30.7 

0.153 

30.2 

0.201 

30.2 

0.211 

30.4 

The  values  of  s  and  c.e.  obtained  when  the  lower  basic 
slope  is  positive  fluctuate  greatly  but  when  the  lower  basic 
slope  becomes  negative  the  c.e.  becomes  almost  constant,  about 
30,  and  the  skewness  stabilizes  near  0,2.   Therefore,  it  is 
concluded  that  the  model  of  a  temperature  signal  has  a  lower 
basic  slope  less  than  zero. 
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-0.6 
FIGURE  16: 


-0.2         0 

LOWER  BASIC  SLOPE 

VARIATION  OF  COMPLEX  RAMP  SIGNAL  WITH  LOWER  BASIC 
SLOPE  (PERIOD  =  800,  UPPER  BASIC  SLOPE  =  0.8, 
RAMP  WIDTH  =  40) 


TABLE  XII. 


VARIATION  OF  DIFFERENTIATED  COMPLEX  RAMP  SIGNAL 
WITH  RAMP  WIDTH  (FIGURE  17) 


PERIOD 


UPPER  BASIC 
SLOPE 


LOWER  BASIC   RAMP       S. 
SLOPE     WIDTH 


C.E. 


800 


0.8 


-0.4 


40 

0.201 

30.2 

50 

0.145 

23.8 

60 

0.111 

19.5 

70 

0.089 

16.3 

80 

0.073 

14.0 

90 

0.061 

12.1 

LOO 

0.052 

10.6 

The  effect  of  varying  the  ramp  width  has  much  the  same 
effect  as  varying  the  period.   Decreasing  the  ramp  width 
causes  s  and  c.e.  to  increase;  likewise  an  increase  in  period 
has  similar  results.   These  two  adjustments  however  do  not 
produce  equal  changes  in  s  and  c.e.   For  example,  in  TABLE  IX 
there  is  a  signal  with  skewness  of  0.145  and  a  c.e.  of  14.3, 
while  from  TABLE  XII  a  signal  with  skewness  of  0.145  has  a 
cae.  of  23 .8. 
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TABLE  XIII. 


VARIATION  OF  DIFFERENTIATED  COMPLEX  RAMP  SIGNAL 
WITH  PERIOD  AND  RAMP  WIDTH 


PERIOD 

UPPER  BASIC 

SLOPE 

200 

400 

800 

0.8 

1600 

3200 

LOWER  BASIC    RAMP 
SLOPE     WIDTH 


S. 


C.E. 


-0.4 


10 

0.702 

28.1 

20 

0.389 

30.3 

40 

0.201 

30.2 

80 

0.099 

28.5 

160 

0.047 

24.8 

TABLE  XIII  is  listed  for  the  information  of  the  reader  and 
as  an  aid  to  those  who  might  want  to  model  temperature  signals. 
In  each  successive  entry  in  TABLE  XIII  the  period  and  ramp 
width  have  been  doubled.   As  previously  noted,  this  has  the 
effect  of  doubling  the  frequency  of  the  internal  teepees. 
Also,  it  is  a  means  of  changing  the  skewness  while  keeping 
the  c.e.  relatively  constant. 

The  preceding  tables  have  not  been  an  attempt  to  list  all 

the  variations  possible  with  the  different  models,  but  to 

examine  how  various  parameters  affect  the  models  and  to 

determine  which  model,  if  any,  is  the  best  representation  of 

atmospheric  temperature  signals.   As  previously  noted,  the 

complex  ramp  signal  was  selected  as  the  best  representation. 

Obviously  there  exists  an  infinite  number  of  possible  variations 

and  given  the  actual  values  of  a  signal  together  with  the 

previous  tables  one  should  be  able  to  model  most  of  the 

atmospheric  temperature  signals  encountered  in  an  unstable 

atmosphere. 

47 


The  complex  ramp  provides  a  coarse  and  a  fine  "tuning" 
mechanism  for  the  model,.   The  ramp  width  and/or  the  period 
can  be  adjusted  to  the  desired  c.c  value  with  an  s  value 
near  that  desired  and  then  the  skewness  tuned  by  adjusting 
the  upper  basic  slope  while  the  coefficient  of  excess  remains 
constant.   Additionally,  the  model  of  the  differentiated 
complex  ramp  (Figure  6)  is  remarkably  similar  to  the  actual 
differentiated  signal  (Figure  1). 
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V.   CONCLUSIONS 

The  results  of  this  study  indicate  successful  accomplish- 
ment of  the  objective;  which  was  to  model  a  temperature  signal 
from  its  statistical  parameters,  skewness,  and  coefficient 
of  excess.   Furthermore,  additional  insight  was  gained  into 
the  actual  structure  of  temperature  signals. 

Although  the  study  was  not  exhaustive  it  does  indicate 
that  one  cannot  attribute  the  observed  values  of  skewness  and 
coefficient  of  excess  in  the  differentiated  temperature  signal 
entirely  to  the  external  dimensions  and  frequency  of  the  ramps. 
Instead  one  must  look  at  the  internal  structure  of  the  ramp 
itself.   The  ramp  is  composed  of  a  great  many  fluctuating 
high  frequency  saw-teeth  which  gradually  increase  in  amplitude 
giving  an  overall  impression  of  a  relatively  gradual  rising 
ramp.   The  fluctuations  are  presumably  due  to  mixing  or 
entrainment  by  the  diffuse  leading  edge  of  the  manifesting 
plume.   The  saw-teeth  actually  extend  below  the  ambient 
temperature  level.   This  is  not  readily  apparent  when  examining 
temperature  strip  charts  but  can  be  attributed  to  the  mechanical 
response  of  the  recorder  and/or  sensing  system.   The  internal 
saw-teeth  cause  the  variation  in  skewness  while  the  width  or 
frequency  of  the  entire  ramp  causes  variation  in  c.e. 
Therefore,  given  the  skewness  and  coefficient  of  excess  of  a 
temperature  signal,  a  model  can  be  constructed  with  the  exact 
statistics  of  the  original  signal. 
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VI.   RECOMMENDATIONS  FOR  FURTHER  STUDY 

A  natural  extension  of  this  thesis  would  be  an  attempt 
to  model  actual  signals.   This  would  probably  require  not 
only  an  examination  of  strip  chart  recordings  but  also  an 
investigation  of  the  signal's  digital  representation. 

While  the  asymmetrical  ramp  signal  proved  unsatisfactory 
as  a  temperature  model,  in  an  extreme  variation  it's 
statistics  were  remarkably  similar  to  those  of  a  velocity 
signal.   When  the  back  slope  width  was  only  slightly  smaller 
than  half  the  entire  ramp  width  a  small  negative  skewness 
and  a  c.e0  of  near  20  was  noted.   Several  of  these  results 
are  listed  in  TABLE  XIII.   Except  for  a  shape  that  more 
resembles  a  teepee  than  a  ramp,  this  signal  should  behave 
exactly  like  the  previously  noted  asymmetrical  ramp  signal. 
An  inspection  of  a  velocity  signal  indicates  that  this  teepee- 
like representation  is  not  a  bad  resemblance.   However,  it 
should  also  be  noted  that  the  same  results  can  be  achieved 
with  a  complex  ramp  signal  whose  lower  basic  slope  is  greater 
than  zero  but  necessarily  less  than  the  upper  basic  slope 
(TABLE  XI).   Perhaps  a  more  accurate  model  would  be  a  combi- 
nation of  the  two,  where  high  frequency  fluctuations  are 
inscribed  between  two  asymmetrical  ramps  similar  to  those 
noted  above  such  that  the  ramps  have  the  same  widths  but 
different  heights,  i.e.,  front  slopes  and  back  slopes. 
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If  this  model  proved  successful  it  would  hopefully 
inherit  the  characteristic  of  the  complex  ramp  whereby  the 
skewness  and  c.e.  could  be  adjusted  independently. 
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TABLE  XIV. 


VARIATION  OF  A  DIFFERENTIATED  VELOCITY 
SIGNAL  WITH  RAMP  WIDTH 


PERIOD 

BACK 
SLOPE 

FRONT 
SLOPE 

RAMP 
WIDTH 

BACK  SLOPE 
WIDTH 

S. 

C.E. 

1.5 

k 

0.0 

22.0 

100 

-L.5 

1.0 
0.75 

5 
6 

2.0 

-1.53 
-2.89 

20.3 
22.0 

0.60 

7 

-3.59 

24.1 
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( SKEW)    AND    COEFFICIENT 
RAMP    SIGNAL    WITH    A    PERIOD    OF 
G. 8t     A    LOWER    BASIC    SLOPE    OF 
PARAMETERS    BOXED    WITH    *»S    AR! 
A    DO    LOOP    INSERTED    BEFORE 
STATEMENT    90    WILL    ALLOW    ONE 
INCREMENTED. 


THIS  PROGRAM  COMPUTES  THE  SKEWNESS 
OF  "XCESS  (CE1  FOP  A  COMPLEX 
800 1  AN  UPPER  BASIC  SLOPE  OF 
-0.4  AND  A  RaMP  WIDTH  Oc  40. 
THF  VARIABLES  OF  TH^  SIGNAL. 
STATEMENT  5  AND  LOODED  THPOUGH 
CR    MORE    OF   THE   VARIABLES    TO    BE 

INITIALIZE    DATA 

REAL''4    LSLOP? 
DIMENSION    X( 2  0000) 
5    DO    10    1=1,10000 
10    X(I)=0.0 
NO=100QO 
M  =  l 
CODE =10** 10 

*    INT=80C  '  * 

*  USL0PE=0.8  * 

*  LSL0PE=-0.4  * 

*  WIDTH=4C  * 

IND=0 

Nl=2 

N2=WIDTH/2 

COMPUTE  COMPLEX  SIGNAL 

DO  50  K=l, 10000 

CCMDIJTE    POSITIVE    SLD^~S 

X(1+IND )=2.*USL0P£ 

CO    2  5    J=N1,N2 
25    X  (  2*'  J-  1  + 1 N  0 )  =2  .  ~  J4  U  S  LO  P*7  -  ( 2  .  "  J  - 1  )  *  LS  L 0  P  E 

CCMPUTE    NEGATIVE    SLOPES 

N3=N2-1 
X(2+IND)=3.-*LSL0P£-2.*USL0PE 

DO    35    J=N1 ,N3 
35    X(  2*J+IND)=(  2.  "J  +  D-LSLOPE-2.^  JAUSLOPE 
X(2*M2  +  IND)  =  -2.vN2  'USLOPE 
IND=IND+INT 

CH^CK    IF    SIGNAL    ARRAY     IS     FILLED 

IF( 2*N2+IND.GT. 10000)     GO    TO    75 
50    CONTINUE 


75 
80 
85 
90 
200 


IER  ) 


CALL    MI  SP(  NO,M,X,CODe.  XBAR  ,  STD,  SKP"W,CE  »R  »  N  ,A  ♦  8  ,S  , 
WRIT5:     (6,200     INT, USLOPE,LSLOPE, WIDTH, SKEW, CE 
r.nNTTM!!^ 

X 

T=  ',  I5,2X,  'USLCPE=  ',  F7.3,  2X,  'LSLOP^=', 
=  '  ,F  5.1  ,2  X,1  SKEW=»  ,F3.3,2X  ,  •  C'-  =  *  ,  =8.3, 
TA  SIGNAL  FOLLOWS' ,/,lCO( 10F10. 3/) ) 


CONTINUE 
WRITE  (6,300)  X 
FORMAT  CO'  ,  '  INT  = 
1F7. 3,2X,' WIDTH 
300  FORMAT  (IX. 'DA 
STOP 
END 


/) 
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